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Abstract. Remarkable anisotropic structures have been recently observed in the order parameter Zik of the 
underdoped superconductor Bi2Sr2CaCu208+i • Such findings are strongly suggestive of deviations from a 
simple d^2_y2-wave picture of high-Tc superconductivity, i.e. Zik '~ coskx — cosfc^. In particular, flatter 
nodes in Zik are observed along the kx = ±fcy directions in k-space, than within this simple model for a 
d-wave gap. We argue that nonlinear corrections in the k-dependence of Z\k near the nodes introduce new 
energy scales, which would lead to deviations in the predicted power-law asymptotic behaviour of several 
measurable quantities, at low or intermediate temperatures. We evaluate such deviations, either analytically 
or numerically, within the interlayer pair-tunneling model, and within yet another phenomenological model 
for a d-wave order parameter. We find that such deviations are expected to be of different sign in the two 
cases. Moreover, the doping dependence of the flatness of the gap near the nodes is also attributable to 
Fermi surface effects, in addition to possible screening effects modifying the in-plane pairing kernel, as 
recently proposed. 

PACS. 74.25.-q General properties; correlations between physical properties in normal and superconduct- 
ing states - 74.25.Jb Electronic structure - 74.20.Mn Nonconventional mechanisms - 74.72.Hs Bi-based 
cuprates 

1 Introduction = 0) and the topological nature of the nodal manifold, 

viz. a collection of points, of line segments, or of surface 



patches 



O- Power laws m the low-temperature asymptotic behaviour ii i • r i- i ^ ii • 

, „ 1, ■ -1 C)n the basis ot group theoretical arguments, the sim- 

,^ 01 several linear response electronic properties provide com- i-i-r , r ,■ t 

. r,; 1 piest choice tor a d-wave gap function on a square lattice 



plementary evidence tor d-wave symmetry ot the order pa- : ^ a l^ \ i a ■ rn j j ^ ^ j 

A ri-im l i rnra n IS Zik = /la(k), where Zi IS a J -dependent parameter, and 

rameter (OP) Zlk of high-Tc superconductors as weU aw; 

as preliminary evidence of 'exotic' shapes in the OP of ^ 



5^ ■ heavy fermion superconductors, such as UPts g . This has (jr(k) = -(cos fc^ — cos /cj^) (1) 
' to be contrasted to an "activated" behaviour oc exp(— /3Z\min), 

appropriate of s-wave superconductors, or, in the case of jg the first basis function associated with the d-wave ir- 

mixed symmetry, of superconductors with a non- vanishing reducible representation of the appropriate crystal point 

s-wave contribution to their OP, where Z\mi„ = miuk |Zik| > group, Cav 0. We remark that g{k) is generated, together 

0. In the case of a non-empty nodal manifold for the super- ^^^h an extended s-wave term proportional to 
conducting excitation spectrum iJk, defined as the locus 

of points in k-space such that = 0, a large number of 1 

quasiparticles can be created near such nodes, thus dom- '''(^) ~ 2^*'*^^''^ ^os/cj^), (2) 
inating all the low-temperature electronic properties 

An exact analysis allows one to relate the exponent of the by a nearest-neighbour interaction term in real space. Here 

leading power of the low-T expansion of a given linear re- and in the following we shall measure the wavevectors in 

sponse function to the dimension of the Fermi manifold units of the appropriate inverse lattice spacings. Propor- 

(defined as the locus of states in k-space with vanishing tionality to Eq. (||) allows Z\k to vanish linearly at a given 

dispersion relative to the Fermi level in the normal state, point along the Fermi line, which for most cuprate super- 
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conductors can be modelled by the tight-binding expan- 
sion: 

<^k = — 2t(cos kx + cos ky) + At' cos k^ cos ky — fi = 0, 

(3) 

where t = 0.25 eV, t' = 0.45i measure nearest and next- 
nearest neighbour hopping, respectively, and /i is the chem- 
ical potential. 

On the other hand, increasing experimental evidence 
above all from angle-resolved photoemission spectroscopy 
(ARPES) suggests a richer structure in k-space for the 
OP of the underdoped high-Tc superconductor Bi2Sr2Ca- 
Cu208+5 |6|. In particular, the superconducting gap near 
the nodal points turns out to be flatter than predicted by 
the simple assumption Z\k oc ^(k) |6|. Such a feature is 
consistent with the observation of whole ungapped seg- 
ments of the Fermi line above Tc in the pseudogap regime 
of underdoped Bi2Sr2CaCu208+5 0, and will of course 
serve as a constraint for a microscopic understanding of 
the pairing mechanism. 

Quite remarkably, qualitatively similar deviations from 
a 5(k)-like dispersion have been evidenced in the k-depend- 
ence of the antiferromagnetic gap in the related insulating 
compounds X2CUO2CI2 (X = Ca, Sr) Such a finding 
has been interpreted in terms of an interrelation between 
the antiferromagnetic phase of the parent insulator and 
the underdoped regime of the intervening superconductor 

In this paper, we argue that such extended structures 
in the superconducting OP, interpolating between point 
and line nodes, can be included in the definition of Z\k as 
higher order terms in g(k). We shall then look for their 
signatures in the low-temperature asymptotic electronic 
properties of the superconducting cuprates, as corrections 
to the predicted power-law behaviour. In deriving our re- 
sults analytically, we will specifically consider the inter- 
layer pair-tunneling (ILT) mechanism of high- Tc super- 
conductivity ijl^ , which has been shown to accurately re- 
produce most of the observed gap features pl| . 



2 Extended d-wave gap within the ILT model 

A distinguishing feature of the ILT mechanism, compared 
to other proposed models of HTSC, is that superconduc- 
tivity is driven by a gain in kinetic, rather than potential, 
energy as temperature is lowered below the critical tem- 
perature Tc- It is assumed that coherent single particle 
hopping between adjacent CUO2 layers in the cuprates 
is suppressed by the non-Fermi liquid character of the 
normal state {e.g. due to spin-charge separation), while 
interlayer coherent tunneling of Cooper pairs is allowed 
as soon as a superconducting condensate is established. 
Confined coherence [|2| within CUO2 layers in the normal 
state is indeed largely motivated by the absence of coher- 
ent transport along the c-axis, whereas a comprehensive 
theoretical understanding of it is still lacking. However, 
there is now abundant experimental evidence that c-axis 
transport in the normal state indeed is incoherent, while 



that in the superconducting state may not be p^ . This 
seems to warrant attention being paid to unconventional 
models of high- Tc superconductivity based on relieving c- 
axis frustrated kinetic energy. Recent findings sug- 
gest, however, that the ILT mechanism alone is not suffi- 
cient to account for the large condensation energy Ec, as 
extracted experimentally from measurements of the pen- 
etration length Ac of several single layered compounds, 
such as Tl2Ba2Cu06+5 whereas the predictions 

of the ILT model agrees with the measured value of Ec 
for La2-;rSr;j,Cu04 ||l^,|l^,^|T^. It has been pointed out, 
however, that while considerable experimental effort has 
been devoted to the determination of Ac, extracting Ec 
from existing data on electronic specific heat is by no 
means straightforward |^ . A direct evaluation of Ec from 
its mean-field expression at T = would relieve the 
complications arising from thermal fluctuation effects, in- 
herent in the method of integrating specific heat data, 
from T = through Tc, recently pointed out in Ref. [po| . 
By utilizing the gap equation, Eq. (^), within the ILT 
model, and of the expression relating Ac at T = to Ec 
[|2), we find results for Ac(T = 0) in Bi2212 which are 
within factors of order unity from the experimental values, 
rather than factors of 10 to 20 ||2^. The observed do ping 
dependence of Ac is also qualitatively reproduced [^Bj ^ 

The emerging scenario suggests therefore that some in- 
plane effective interaction might co-operate with the ILT 
mechanism in establishing the superconducting state [p^ . 
One could think of such a mechanism as a seed for the 
Cooper instability, and the origin of the gap's dominant 
d-wave symmetry. Once Cooper pairs are formed in the 
appropriate symmetry channel(s) via such in-plane effec- 
tive interaction, the ILT mechanism would allow the con- 
densate for an additional energy gain, by releasing the 
constraint of in-plane segregation. 

Without explicitly specifying the microscopic origin 
of the in-plane mechanism, we therefore assume the in- 
plane pairing potential to be given by Hck' = ^ff(k)(7(k') 
{V < 0), thus allowing for d-wave symmetry of the order 
parameter. The issue of the competition with other sub- 
dominant (s-wave) symmetry channels in the presence of 
ILT has been addressed in Ref. showing that the d- 
wave contribution wins out at optimal doping and in the 
underdoped regime. Despite its kinetic nature, ILT can 
be absorbed in the interacting part of the Hamiltonian 
as an effective term Tj (k)(5kk', whose k-space locality en- 
forces in-plane momentum conservation during a tunnel- 
ing process 0. Following Ref. |l^, we assume Tj(k) = 
t']_(k)/t = Tjg^(k), being t±{k) the single-particle inter- 
layer hopping amplitude, with i±(k) cx 5^(k), as suggested 
by ARPES as well as by band structure calculations 
p5[ . A standard mean-field diagonalization technique then 
yields the following expression for the energy gap p6|,p|: 



l-Tj(k)xk 



(4) 



where Xk = (2ii'k) ^ tanh(/3i?k/2) is the superconducting 
pair susceptibifity, and T'k = (Ck + l^kP)^^^ is the up- 
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per branch of the superconducting elementary excitation 
spectrum. 

Along the Fermi line (^k = 0) at T = 0, one immedi- 
ately sees that: 

Z\k = Ag{k) + ^Tjg^ik) sgn[g{k)]. (5) 

Such an expression, together with manifestly fulfilling the 
requirement of d-wave symmetry, also endows the super- 
conducting gap with a richer structure near the nodal 
points along the — iky directions. This is probably 
best seen by considering the Fourier expansions: 

oo 

5(k) = -2 ^ J4m-2{k) cos[(4m - 2)4^], (6a) 

m — 1 

oo 

h(k) = Jo(fc) + 2 ^ Jimik) cos(4m0), (6b) 

m— 1 

Tjik)/Tj = — + ^+ ^ a4™cos(4m0), (6c) 

m— 1 

with 

a4m = ^^4m(4A;) + ij4,n(2fc) + ^ l)"->^4m (2fc\/2) 

62 2 lb 

-b-iy^J^^{kV2) - ]ji,n (^^) cos(4m0o). 
4 4 \ sm 00 / 

(7) 

Here, the generic wavevector k is expressed in terms of its 
modulus k and of the angle cj) formed with the FX direc- 
tion in the first Brillouin zone (IBZ), k = (fc cos 0, k sin (/)), 
Ja{x) are Bessel functions of the first kind and order a, 
and tan 0p = 5 ■ 

Eq. (H) is to be contrasted to the phenomenological fit 

Z\k = A[B cos(20) + {l-B) cos(6(/))] (8) 

proposed in Ref. for Z\k along the Fermi line: Instead 
of requiring an in-plane interaction extended to further 
neighbours, Eq. (||) endows the superconducting gap with 
the observed flat structure around the nodes, through the 
ILT term Tj{k). In Fig. 0, we fit Eq. (|) against Mesot et 
al.^s experimental data for one of the underdoped Bi2212 
samples in Ref. having Tc — 75K. A remarkable agree- 
ment follows already by fixing A so that |Z\k| reproduces 
the maximum datum at k = (tt, tt), whereas Tj is taken 
to be 0.04 eV jlO). In particular, besides obtaining an en- 
hanced maximum value of |Z\k| at k = (tt, tt), we are thus 
able to recover the anomalously flat region around the 
node at 4> = 45° in a rather natural way. We note, how- 
ever, that our fit requires Z\ « T7/2 around optimal dop- 
ing, which will not be without consequences in evaluating 
other fundamental quantities [23| . 

Eq. (||) already contains the doping dependence of the 
observed gap anisotropy, although in a hidden way. As 
pointed out in Ref. the auxiliary parameter A is to 
be self-consistently determined by solving the appropriate 
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Fig. 1. Fit for |Zik| within the ILT model, Eq. (||) (solid 
line), and in the case of a simple d-wave gap |zik| = ^5(k) 
(dashed line), against Mesot et al.'s ARPES data for under- 
doped Bi2212 {Tc = 75 K, Ref. §). 

gap equation. Besides being intrinsically doping depen- 
dent, this equation is unconventionally modified by the 
presence of a k-local effective interaction, as induced by 
the ILT mechanism. Moreover, the role of the contribution 
cx Tj(jr^(k) in Eq. (||) is strongly influenced by the actual 
location of the Fermi line, as g*{k) is sharply peaked at 
k = (OjTf) (and symmetry related points). 

Eq. (^ also facilitates the evaluation of the slope of 
the superconducting gap — {l/2)d\Ak\/d(j) at the nodal 
point along the Fermi line. Such a quantity is related to the 
temperature derivative of the superfluid stiffness at T = 0. 
In particular, it is seen that the ratio u/i/Zimax decreases 
with underdoping From Eq. (|^), one derives that va 
is independent of Tj, and that therefore a doping induced 
change of w/i through A essentially can be traced back to 
the actual position of the Fermi line, as discussed above, 
within the ILT model. The ratio W/i/^max will anyway 
deviate from its value within simple d-wave (BCS-like) 
models, as a function of doping, due to the enhancement 
of zimax induced by ILT. 

3 Low-temperature asymptotic behaviour of 
electronic properties 

We now address the issue, whether such extended fea- 
tures of the OP near the nodes, as those described in the 
previous section, induce deviations in the low or inter- 
mediate temperature asymptotic behaviour of linear re- 
sponse electronic properties in the superconducting state. 
In what follows, we shall limit our discussion to clean 
superconductors, and neglect impurity effects altogether. 
Mean-field (BCS or BCS-like) expressions for most lin- 
ear response electronic properties are available also in the 
case of anisotropic, i.e. non s-wave, superconductors. In 
particular, we have in mind observable quantities such 
as the superconducting density [0, the electronic spe- 
cific heat the spin susceptibility the penetration 
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depth 1 28 , the thermal conductivity ||2^ , and so on. Their 
expressions basically involve the evaluation of some inte- 
gral of the kind: 



(9) 



where (3 = {k^T)"^, fk{P) is a (dimensional) function of 
wavevector k and temperature, related to the electronic 
quantity of interest, and the integration is extended to 
the IBZ, k € [— TT, tt] x [— tt, tt] (see App. |^. In the case 
of d-wave superconductors, i?k is allowed to vanish at the 
intersection between the Fermi line and the nodal lines of 
the gap function. Around such points, quasiparticles can 
be created in large numbers. In the limit of low temper- 
atures {P — *■ oo), therefore, the value of the integral in 
Eq. (|^) is dominated by the contributions from wavevec- 
tors k close to such point nodes. Around such nodes, it is 
useful to introduce the new sets of coordinates (fci, or 
{e,0), defined as [|: 

— vp • k = vpki = ecosd, (10a) 
Z\g(k) ~ V2 • k = ■(72^2 = e sin 61, (10b) 

in units where h = I. Here, vp and V2 are the Fermi ve- 
locity and a suitable 'gap' velocity, respectively, evaluated 
at i?k = 0, and e measures the distance in energy from a 
given dispersionless point implicitly defined by £'k — 0. In 
terms of the new coordinates, the superconducting spec- 
trum for a simple d-wave superconductor near a node 
therefore looks like an anisotropic Dirac cone W, 



v^kl 



V2 k2 



1/2 



(11) 



The observation of flatter structures near the nodes [g| 
not only implies a more significant anisotropy ratio vy/v2, 
but also the possibility that higher order terms in e may 
contribute to i?k, Eq. (11). Indeed, within the ILT model, 
from Eq. (0) at T = one obtains 



(12) 



to lowest order in e/e^, with l/el {l/2){Tj / A^) related 
to the pair-tunneling arnplitude Tj and to the auxiliary 
gap parameter A (Figs. Q and ^). 

Other models, based on extended in-plane pairing mech- 
anisms, would in general yield different polynomial cor- 
rections in e to E\^. For instance, within the spin fluctua- 
tion theory , the following phenomcnological expansion 
holds for the momentum distribution of the superconduct- 
ing energy gap Q 



N 



zik = ^5(k)^d„/i"(k). 



(13) 



n=0 



with all coefficients dn = 1. We explicitly observe that 
for N ~ 0, do ~ I, one recovers the simple d-wave gap 




Fig. 2. Typical contour lines of the superconducting spectrum 
Ek in the simple d-wave case (dashed lines) and in presence of 
ILT (continuous line). 



Z\k oc g(k), while the case = 1, with the identifications 
A ^ BA, do = l,di= 4(1 - B)/B, maps to i,|l|: 



Au ^ A[Bg(k) + (1 - B)g{2k)], 



(14) 



which is compatible with the phenomcnological fit Eq. (|^) 
proposed by Mesot et al. in Ref. for their experimen- 
tal data of |Z\k| along the Fermi line 1^,0. In partic- 
ular, Eq. (^_4|) would follow from a correction (5Vkk' oc 
.g(2k)(7(2k') to the in-plane coupling, corresponding to 
next-nearest neighbours interaction. 

In such a particular case, and assuming for simplicity 
t' ~ fi — in Eq. (||), one straightforwadly obtains 



Eu 



cos^ + sin^ \ 1 cos ( 



1/2 



(15) 



where = tB / {1 ~ B) is now related to the ratio of near- 
est vs next-nearest neighbours coupling. Therefore, both 
within the ILT model and within other models, based on 
extended in-plane pairing, the additional mechanism re- 
sponsible for the nonlinear correction to E]^ away from its 
nodes introduces new energy scales (here, or e*, respec- 
tively). Fig. 1^ depicts the two different ways in which E]^ 
deviates from the cone-like shape, Eq. (|ll|), near a node, 
in the two cases given by Eqs. (|l2) and ([l5[). 

In the absence of any such additional mechanism (e^,, e^, = 
0), the leading contribution to Eq. (^) for the simplest, 
reference case (/3k(/3) = 1 is: 



(16) 



where A = {2ttvpV2A)~^ is a doping-dependent factor, 
and = denotes equality up to terms vanishing exponen- 
tially with P at all energy scales, as /3 00 (T ^ 0). 
Eq. (^6|) should be regarded as typical of the power-law 
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Fig. 3. Deviations from the simple d-wave case, Eq. (^l|) (con- 
tinuous line), of the superconducting spectrum around a 
node, within the ILT model, Eq. (p^ (dashed line), and in the 
case of an extended d-wave gap, Eq. (|]_g) (dashed-dotted line), 
as a function of the reduced coordinates e/e^ (^Z?*) respec- 
tively). Such deviations are most easily seen along the direc- 
tion of the nodal line (left panel, V2k2 — or 6 — n), and along 
the Fermi line (right panel, wpfci — or 8 = -k /2). Note that 
along V2k2 = 0, one has ~ e also within the ILT model. 



asymptotic low-temperature behaviour of the supercon- 
ducting electronic properties within a simple d-wave BCS- 
like model. 

In order to obtain an asymptotic expansion for !Fi{j3) 
as /3 ^ oo (T ^ 0), including the corrections due to ILT, 
Eq. (^2|), we observe that the integration over e in Eq. (||) 
is actually made of two contributions: 



dt 



de 



de. 



(17) 



In the first integral, we may safely retain only the linear 
term ~ e in the exponent, since e < e*. In the second 
contribution, this is no longer possible, and Eq. (12) has 
to be retained in full. However, since e > > 0, one 
can make use of Laplace's (saddle point) method for the 
integral over angles around 9 = 0. The final result is: 



A 
A 



1.4 
1.2 
1 



simple d-wave 
ILT 

extended d-wave 




0.4 0.6 

Fig. 4. Asymptotic power-law (oc T^, solid line) and modified 
power-law behaviours of J-'i (/3) , in the presence of ILT (dashed 
line), and in the extended d-wave case (dashed-dotted line), as 
(/3e*)~^ — > 0. Given the values used for the fits of |zik| along 
the Fermi line in Refs. ^ and it turns out that ~ 250 K. 



{2TrvpV2BA)^^ and e^, i— > e^. In this case 
vides i?k with a different kind of anisotro 
to the simple case, Eq. (pj]), than Eq. (|l 
in the latter case one always has > e. 



1 - (1 + f3e,)e''''* + ( ^, /?£ 



1 



1- 1 



3 



where r{x), r{a,x) are Euler gamma and incomplete 
gamma functions, respectively [p2|[. A comparison of Eq. ( p^ 
and Eq. (jl^) is provided by FigT^ and shows that !Fi{p) 
gets effectively suppressed in the presence of flat nodes in 
the order parameter, as provided by the ILT mechanism, 
with respect to the simple d-wave case, at an energy scale 

No such simple asymptotic expansion for J-i (/3) is avail- 
able in the extended d-wave case described by Eq. (p^, 
and the integrations have to be performed numerically. 
Fig. U shows the result, with the identifications A^ A = 



Eq. (|T|) pro- 
with respect 
does. While 
here one has 

i?k ^ e, depending on the angle 6 (cfr. Fig. |3|). As a con- 
sequence, Ti{l3) is enhanced with respect to the simple 
d-wave case, at an energy scale ~ e*. 



4 Conclusions 

Motivated by recent experimental findings of extended flat 
structures in the order parameters of the underdoped d- 
wave superconductor Bi2Sr2CaCu208+5 1^, we have ad- 
dressed the issue of whether nonlinear, high-energy correc- 
tions to the superconducting energy spectrum around 
the gap nodes induce deviations in the predicted power- 
law behaviour of several electronic properties at low or in- 
termediate temperatures. We have shown that nonlinear 
corrections to E\^ in general introduce additional energy 
scales in the problem. Deviations from the usual power- 
w behaviour of the superconducting electronic proper- 
i ;s are indeed to be expected at such energy scales, but 
the actual value and sign of such deviations are specific 
(18) the model under consideration. In particular, within 
the ILT model, we have explicitly derived the expected 
corrections to a typical power-law asymptotic behaviour 
as T — > 0, showing these to be negative, whereas within 
a phenomenological model of extended d-wave supercon- 
ductivity 1^ such corrections are predicted to be positive. 
Whether such deviations will actually be observable in 
real measurements of superconducting electronic proper- 
ties, will of course depend on the effective values of the 
additional energy scales e^, or e^, in real compounds. 
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A Low-temperature superconducting 
electronic properties 

We now give a sketch of how the low-temperature asymp- 
totic behaviour of several electronic properties of interest 
can be reduced to that of Ti{f3) or its derivatives. Most 
electronic quantities in the superconducting state are in 
fact given by Eq. with (/?k(/3) actually depending on k 
only through E\^. In what follows, /(e) = [1 + exp(/3e)]~^ 
denotes the Fermi function. 

Within BCS theory, the electronic specific heat is given 
by H: 



Cv =^2fcB/3£^k 



dfi 



Of \ 
dEi,) 



k 



2kBf3E^ - 



dEu 



2fcB/3'^r(/3), 



(19) 



where apices denote derivatives with respect to f3. Here, 
</'k(/3) = 2kBf3^E^, and we have made use of the fact that 
{-df/dE^)=l3exp{(3Eu). 

Analogously, the unrenormalized, static, isotropic spin 
susceptibility xo = Xo(q ^ 0,ti; ^ 0), which is directly 
related to the Knight shift, is simply given by |pT|,^: 



Xo = 



E 



9^k 



The expression of the electronic thermal conductivity 
for an anisotropic d-wave superconductor also involves an 
average of {-df/dEi,) over the IBZ 

2 



dE. 



Ei 



OEy 



r(k), (21) 



where T(k) is the superconducting quasiparticles lifetime. 
Due to the presence of the x component of the group veloc- 
ity Vk-Ek, however, its expression in our notation reduces 
to: 

_ J_4 j_ f°° . f^"" 
87r2 V2 T Jo Jq 



det' I d^i-^j ^cose*-!- ^sine* 



B A limiting case 

In the absence of in-plane coupling, a spurious solution of 
the mean-field gap equation at T can be implicitly 
expressed via [nil : 



Eu 



-T,(k) = -Tjg\k). 



(24) 



In such a limiting case, the superconducting energy spec- 
trum would have purely kinetic origin, and would be iden- 
tified with the interlayer pair-tunneling amplitude, divided 
by two. A closed expression can then be obtained for 
J-'i{(3), by utilizing the useful result: 



1 



(2^)2 



dxg{x)K{\/l- x^), (25) 



where t/[?7(a;)] is any continuous functional of //(k) = /i(k) 
or g(k) alone, and K{x') [x' = ^/l — x'^) is the complete 
elliptic integral of first kind |^ . From Eq. (^) , expanding 
K{x') around x = one eventually arrives at the closed 
expression: 



1 



4 ; 



21og2-iv Q) + logC'/^ 



(26) 



where ip{x) is the digamma function js^, and the ILT 
amplitude Tj itself here fixes the appropriate energy scale, 
through C = \[3Tj. 



(20) References 



(22) 

where £o = i'FT(kF) is the quasiparticle mean free path 
at the nodes. The final result crucially depends on the 
anisotropy ratio V2/vp, and would be different in the two 
cases given by Eqs. (^2|) and (^5|), due to their different 
9 dependence. This has to be contrasted with the result 
obtained in the simple d-wave case, where |34[1: 



with r] 



Ke = — 1 + ^ 

{Stt)-^ dxx^-df/dx). 



(23) 
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